1 


Climate Hazard Assessment for Stakeholder 


2 Adaptation Planning In New York City 


3 


4 

5 Radley M. Horton*, Vivien Gomitz*, Daniel A. Bader*, Alex C. Ruane^’*, Richard Goldberg*, 

6 and Cynthia Rosenzweig^’* 

7 


8 

9 * Center for Climate Systems Research, Earth Institute, Columbia University, New York, NY 

10 

11 NASA Goddard Institute for Space Studies, New York, NY 

12 


13 

14 

15 Corresponding author address: Radley Horton, Columbia University Center for Climate 

16 Systems Research, 2880 Broadway, New York, NY, 10025. 

17 Email: rhl42@columbia.edu 

18 


1 



19 


Abstract 


20 This paper describes a time-sensitive approach to climate change projections, developed 

21 as part of New York City’s climate change adaptation process, that has provided decision 

22 support to stakeholders from 40 agencies, regional planning associations, and private companies. 

23 The approach optimizes production of projections given constraints faced by decision makers as 

24 they incorporate climate change into long-term planning and policy. New York City 

25 stakeholders, who are well-versed in risk management, helped pre-select the climate variables 

26 most likely to impact urban infrastructure, and requested a projection range rather than a single 

27 ‘most likely’ outcome. The climate projections approach is transferable to other regions and 

28 consistent with broader efforts to provide climate services, including impact, vulnerability, and 

29 adaptation information. 

30 The approach uses 16 Global Climate Models (GCMs) and three emissions scenarios to 

31 calculate monthly change factors based on 30-year average future time slices relative to a 30- 

32 year model baseline. Projecting these model mean changes onto observed station data for New 

33 York City yields dramatic changes in the frequency of extreme events such as coastal flooding 

34 and dangerous heat events. Based on these methods, the current 1-in-lO year coastal flood is 

35 projected to occur more than once every 3 years by the end of the century, and heat events are 

36 projected to approximately triple in frequency. These frequency changes are of sufficient 

37 magnitude to merit consideration in long-term adaptation planning, even though the precise 

38 changes in extreme event frequency are highly uncertain. 

39 
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40 


1. Introduction 


41 This paper describes a methodological approach to stakeholder-driven climate hazard 

42 assessment developed for the New York Metropolitan Region (Fig. 1). The methods were 

43 developed in support of the New York City Panel on Climate Change (NPCC; NPCC 2010). 

44 The NPCC is an advisory body to New York City’s Climate Change Adaptation Task Force 

45 (CCATF), formed by Mayor Michael Bloomberg in 2008 and overseen by the Mayor’s Office of 

46 Long Term Planning and Sustainability. As described in NPCC (2010), the CCATF is 

47 comprised of stakeholders from 40 city and state agencies, authorities, regional planning 

48 associations, and private companies, divided into four infrastructure workgroups 

49 (communication, energy, transportation, and water and waste), and one policy workgroup. 

50 The CCATF effort was motivated by the fact that New York City’s population and 

51 critical infrastructure are exposed to a range of climate hazards, with coastal flooding associated 

52 with storms and sea level rise the most obvious threat. Approximately 7 % (11%) of NYC area 

53 is within 1 meter (2) of sea level (Weiss et al. 2011). A recent study ranked NYC 7* globally 

54 among port cities in exposed population and 2”‘^ globally in assets exposed to storm surge 

55 flooding and high winds (Nicholls et al. 2008). Furthermore, because NYC, like much of the 

56 U.S. (ASCE, 2009), has aging infrastructure, climate vulnerability may be enhanced. By 

57 showing leadership in the infrastructure adaptation process, the NYC effort may be able to 

58 provide lessons to other cities as they plan adaptation strategies. 

59 Stakeholder input regarding climate information was collected in several ways. Between 

60 September of 2008 and September of 2009, each CCATF sector working group held monthly 

61 meetings in conjunction with the Mayor’s Office of Long Term Planning and Sustainability. 

62 During the initial meetings, representatives from each sector identified key climate hazards; they 
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63 also interacted iteratively with the scientists, seeking clarification, and requesting additional 

64 information. They commented on draft documents describing the region’s climate hazards, and 

65 climate seminars were held with individual agencies as requested. The climate hazard 

66 assessment process was facilitated by prior collaborative experience between the NPCC’s 

67 climate scientists and stakeholders in earlier assessments, including the Metro East Coast Study 

68 (MEC; MEC 2001), as well as work with the New York City Department of Environmental 

69 Protection (NYCDEP; NYCDEP 2008; Rosenzweig et al. 2007) and the Metropolitan Transit 

70 Authority (MTA; MTA 2007). 

71 The climate hazard approach is tailored towards impact assessment; it takes into 

72 consideration the resource and time constraints faced by decision makers as they incorporate 

73 climate change into their long-term planning. Eor example, the formal write-up of the climate 

74 risk information was needed within less than 8 months of the NPCC’s launch (NRC 2009); given 

75 this time frame and the broad array of stakeholders in the CCATE, a standardized set of climate 

76 variables of broad interest were emphasized, with the understanding that future studies could 

77 provide climate information tailored to more unique applications^ 

78 Within this framework, the NPCC worked with stakeholders to pre-select for analysis 

79 climate variables and metrics most likely to impact existing assets, planned investments, and 

80 operations (Horton and Rosenzweig 2010). Eor example, the number of days below freezing was 

81 identified as an important metric for many sectors, due to the impacts of freeze-thaw cycles on 

82 critical infrastructure (New York City Climate Change Adaptation Task Eorce, 2008-2009). Due 

* A tailored assessment of changes in snow depth and timing of snow melt in the Catskill Mountains approximately 
100 miles north of New York City (NYC DEP, 2008) would be of interest to managers of only a small but important 
subset of infrastructure — reservoirs and water tunnels. Such a fme-scale assessment would benefit from more 
complex downscaling approaches than those applied here. 
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83 to the diversity of ageneies, projections were requested for multiple time periods spanning the 

84 entire 21®' century. 

85 Stakeholders also helped determine the presentation of climate hazard information. For 

86 example, because NYC stakeholders are used to making long-term decisions under uncertainty 

87 associated with projections of future revenues, expenditures, and population trends, for example, 

88 they preferred projection ranges to a single ‘most likely’ value (New York City Climate Change 

89 Adaptation Task Force, 2008-2009). 

90 Itemized risks associated with each climate variable were ultimately mapped to specific 

91 adaptation strategies. For example, more frequent and intense coastal flooding due to higher 

92 mean sea level was linked to increased seawater flow into New York City’s gravity-fed and low- 

93 lying Wastewater Pollution Control Plants (WPCP), resulting in reduced ability to discharge 

94 treated effluent (NPCC 2010; NYCDEP 2008). NYCDEP is reducing the risk at the Ear 

95 Rockaway Wastewater Treatment Plant by raising pumps and electrical equipment to 14ft above 

96 sea level based on the projections described here (NYC Office of the Mayor, 2009). 

97 Climate hazard assessment was only one component of the NPCC’s impact and 

98 adaptation assessment. Vulnerability of infrastructure (and the populations that rely on it) to 

99 climate impacts can be driven as much by its state of repair (and how it is used) as by climate 

100 hazards (NRC, 2009). Climate adaptation strategies should be based on many non-climate related 

101 factors, such as co-benefits (e.g., some infrastructure investments that reduce climate risks will 

102 also yield more efficient and resilient infrastructure in the face of non-climate hazards; NRC 

103 2010a) and co-costs (e.g., adapting by using more air conditioning increases greenhouse gas 

104 emissions). NPCC experts in the risk management, insurance, and legal fields provided guidance 

105 on these broader issues of vulnerability and adaptation, developing for example an eight step 
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106 adaptation assessment proeess and templates for ranking relative risk and prioritizing adaptation 

107 strategies (NPCC, 2010). This paper foeuses on the provision of stakeholder-relevant elimate 

108 information in support of the broader NPCC assessment. 

109 Seetion 2 deseribes the methodology used for the NPCC’s elimate hazard assessment. 

110 Seetion 3 eompares elimate model hindeasts to observational results for the New York 

111 Metropolitan Region. Hindeast results are a reeurring stakeholder request, and they helped 

112 inform the global elimate model (GCM)-based projeetion methods. Seetion 4 doeuments the 

113 regional projeetions, in the eontext of stakeholder usability. Seetion 5 eovers eonelusions and 

1 14 reeommendations for future work. 

115 

116 2. Methodology 

117 a. Observations 

118 Observed data are from two sources. Central Park station data from the National Oceanic 

119 and Atmospheric Administration, National Climatic Data Center, United States Historical 

120 Climatology Network (NOAA NCDC USHCN) Version 1 data set (Karl et al. 1990; Easterling 

121 et al. 1999; Williams et al. 2005) formed the basis of the historical analysis and projections of 

122 temperature and precipitation. Gridded output corresponding to New York City from the 

123 National Centers for Environmental Prediction / Department of Energy (NCEP/DOE) Reanalysis 

124 2 output (Kanamitsu et al. 2002) is also used for GCM temperature validation (section 3). 

125 

126 b. Climate Projections: General Approach 

127 1 ) Global C limate Models and emissions scenarios 
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128 Climate projections are based on the coupled GCMs used for the Intergovernmental Panel 

129 on Climate Change Fourth Assessment Report (IPCC; IPCC, 2007). The outputs are provided by 

130 the World Climate Research Programme's (WCRP's) Coupled Model Intercomparison Project 

131 phase 3 (CM1P3) multi-model dataset (Meehl et al. 2007a). Out of 23 GCM configurations from 

132 16 centers, the 16 GCMs that had available output for all three emissions scenarios archived by 

133 WCRP were selected (Table 1; Al, AIB, and Bl; IPCC Special Report on Emissions Scenarios; 

134 SRES; Nakicenovic et al. 2000). 

135 The 16 GCMs and three emissions scenarios combine to produce 48 output sets. The 48 

136 members yield a model and scenario-based distribution function based on equal weighting of 

137 each GCM and emissions scenario. The model-based results should not be mistaken for a 

138 statistical probability distribution (Brekke et al. 2008) for reasons including the following: 1) no 

139 probabilities are assigned by the IPCC to the emissions scenarios^ 2) GCMs are not completely 

140 independent; with many sharing portions of their code and a couple differing principally in 

141 resolution only, and 3) the GCMs and emissions scenarios do not sample all possible outcomes, 

142 which include the possibility of large positive ice-albedo and carbon cycle feedbacks, in addition 

143 to uncertain aerosol effects. Caveats notwithstanding, the model-based approach has the 

144 advantage (relative to projections based on single numbers) of providing stakeholders with a 

145 range of possible outcomes associated with uncertainties in future greenhouse gas concentrations 

146 (and other radiatively important agents and climate sensitivity (NRC, 2010b). 

147 Some authors (see e.g.. Smith et al. 2009; Tebaldi et al. 2005; Greene et al. 2006; Brekke 


^ It has been argued that since high global anthropogenic CO 2 emissions growth rates (3.4 % year * between 2000 
and 2008; Le Quere et al. 2009) led to 2008 estimated emissions reaching the levels of the highest SRES scenario 
(AIFI), other SRES scenarios may be unrealistically low. 
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148 et.al. 2008; Georgi and Mearns 2002) have explored alternate approaehes that weight GCMs 

149 based on eriteria ineluding hindeasts of regional elimate or key physieal proeesses. That more 

150 eomplex approaeh is esehewed here in favor of equal GCM weighting for several reasons. First, 

151 beeause model ‘sueeess’ is often region and variable-speeifie, and stakeholders differ in their 

152 elimate variables and geographieal ranges of interest^, produetion of eonsistent seenarios based 

153 on model weighting is a major researeh effort beyond the seope of New York City’s initial 

154 assessment. Seeond, while long-term researeh eould be geared towards developing optimized 

155 multivariate (and/or multi-region) weighting, researeh suggests that eompensating biases tend to 

156 yield eomparable model performanee (Brekke et al. 2008). Third, historieal aeeuraey may have 

157 been aehieved for the ‘wrong’ reasons (Brekke et al. 2008) and GCM hindeasts did not share 

158 identieal foreing, espeeially with respeet to aerosols (Rind et al. 2009). Fourth, shifting elimate 

159 proeesses with elimate ehange may favor different models in the future. Finally, eliminating 

160 ensemble members reduees the representation of uneertainty relating to elimate sensitivity. 

161 

162 2) Timeslices 

163 Beeause eurrent-generation GCMs used for elimate ehange applieations have freely 

164 evolving oeean and atmospherie states, they are most appropriate for deteetion of long-term 

165 elimate and elimate ehange signals. The 30-year timesliee applied here is a standard timeseale 

166 (WMO 1989) that represents a middle ground, allowing partial eaneellation of eurrently 

167 unpredietable interannual to interdeeadal variability (maximized by ineluding many years), while 

168 maintaining relatively monotonie anthropogenieally-indueed foreing trends (maximized by 

169 ineluding few years). The ‘1980s’ timesliee represents baseline eonditions between 1970-1999; 

^ New York City’s task force included corporations with national and international operations. 
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170 future timeslices for the 2020s, 2050s, and 2080s are similarly defined. 

171 

172 3) C LIMATE CHANGE FACTORS AND THE DELTA METHOD 

173 Mean temperature ehange projeetions are expressed as differences between each model’s 

174 future timeslice simulation and its baseline simulation; mean precipitation is based on the ratio of 

175 a given model’s future to its baseline values. This approach offsets a large source of model bias: 

176 poor GCM simulation of local baseline conditions (section 3b) arising from a range of factors 

177 including the large difference in spatial resolution between GCM gridboxes and station data. 

178 Because monthly averages from GCMs are generally more reliable than daily output 

179 (Grotch and MacCracken 1991), monthly mean GCM changes were projected onto observed 

180 1971-2000 daily Central Park data for the calculation of extreme events"^. This simple and low- 

181 cost downscaling approach is known as the delta method (Gleick 1986; Arne 11 1996; Wilby et al. 

182 2004). Like more complex statistical downscaling techniques (e.g., Wigley et al. 1990), the delta 

183 method is based on stationarity (see e.g., Wilby et al. 1998 and 2002; Wood 2004), and largely 

184 excludes the possibility of large variance changes through time, although for the Northeast U.S. 

1 85 such changes are uncertain^ . 

186 More complex statistical approaches, such as those that empirically link large-scale 

187 predictors from a GCM to local predictands (see e.g., Bardossy and Plate 1992) may yield more 

188 nuanced downscaled projections than the delta method. These projections are not necessarily 

For coastal flooding and drought, the 20* century was used as a baseline, due to high interannual/multidecadal 
variability and policy-relevance of 1 -in- 100 year events. 

^An exception may be short-term precipitation variance, which is expected to increase regionally with the more 
intense precipitation events associated with a moister atmosphere (see e.g., Emori and Brown 2005; Cubasch et al. 
2001;Meehlet al. 2005) 
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189 more realistic, however. Historical relationships between large-scale predictors and more 

190 impacts-relevant local predictands may not be valid in a changing climate (Wilby et al. 2004). 

191 GCM development and evaluation has also historically been more focused on seasonal and 

192 annual climatologies than the daily and interannual distributions that drive analogue approaches. 

193 Table 2 provides a set of stakeholder questions to inform the choice of downscaling technique, a 

194 topic that is discussed further in Section 5. 

195 

1 96 4) Spatial extent 

197 The projections are for the land-based GCM gridbox covering New York City. As shown 

198 in Fig. 2, the 30-year averaged mean climate changes are largely invariant at sub-regional scales; 

199 the single grid box approach produces nearly identical results to more complex methods that 

200 require extraction of data from multiple gridboxes and weighted spatial interpolation. As shown 

201 in section 4e, for the metrics evaluated in this study, the GCM gridbox results also produce 

202 comparable results to finer resolution statistically and dynamically downscaled products. Since 

203 baseline climate (as opposed to projected climate change) does differ dramatically over small 

204 spatial scales (due to factors such as elevation and surface characteristics), and these fine-scale 

205 spatial variations by definition cannot be captured by course-resolution GCMs, GCM changes 

206 are trained onto observed Central Park data using the procedures described above. 

207 

208 5) N UMBER OF SIMULATIONS 

209 For 13 of the 16 GCMs’ Climate of the 20* Century and future AIB experiments, and 7 

210 of the 16 B1 and A2 future experiments, multiple simulations driven by different initial 

211 conditions were available. Analysis of hindcasts (Table 3a) and projections (Table 3b) from the 
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available NCAR CCSM coupled GCM simulations^ revealed only minor variations in 30-year 
averages, suggesting that one simulation per model is sufficient. Using an ensemble for each 
GCM based on all the available simulations with that GCM is an alternative approach; however, 
the effort and data storage needs may not be justified given the similarity of the ensemble and 
individual simulation results shown in Table 3. Furthermore, ensemble averaging unrealistically 
shrinks the temporal standard deviation^. 

c. Climate Projections: Sea Level Rise 

To address large uncertainties associated with future melting of ice sheets, two sea level 
rise projection methods were developed; these are referred to as the IPCC -based and rapid ice 
melt scenarios respectively. 

1) IPCC AR4-BASED APPROACH 

The IPCC Fourth Assessment (AR4) approach (Meehl et al. 2007b) was regionalized for 
New York City utilizing four factors that contribute to sea level rise; global thermal expansion, 
local water surface elevation, local land uplift/subsidence, and global meltwater^. Thermal 
expansion and local water surface elevation terms are derived from the GCMs (outputs courtesy 
of WCRP and Dr. Jonathan Gregory, personal communication). Local land subsidence is derived 
from Peltier (2001) and Peltier’s ICE-5Gvl.2 ice model (2007) 
(http;//www.pol. ac.uk/psmsl/peltier/index.html) . The meltwater term was calculated using mass 
balance temperature sensitivity coefficients for the different ice masses, based on observed 

® This GCM was selected because it provided the most 20* *'’ and 2L' century simulations 

^ This is a general criticism; for the particular case when the delta method is used (as here) shrinking of the temporal 
standard deviation has no bearing on the results 

* Only seven GCMs provided outputs for sea level rise projections; see Horton and Rosenzweig (2010) for 
additional information. 
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233 historic relationships between global mean surfaee air temperature, iee mass, and rates of sea 

234 level rise (Meehl et al. 2007b)^. Regionalization of sea level rise projections, based on the four- 

235 eomponents deseribed above, have been used in other studies (e.g.. Mote et al. 2008). 

236 

237 2) Rapid ice melt scenario 

238 Beeause of large uneertainties in dynamieal iee sheet melting (Hansen et al. 2007; Horton 

239 et al. 2008), and recent observations that iee sheet melting has aeeelerated within this past deeade 

240 (e.g., Chen et al. 2009), an alternative sea level rise seenario was developed. This upper bound 

241 sea level rise scenario allowing for rapid iee melt was developed based on paleo-sea level 

242 analogues, in particular the ~10,000-12,000-year period of rapid sea level rise following the end 

243 of the last iee age (Peltier and Fairbanks 2006; Fairbanks 1989). While the analogue approach 

244 has limitations (most notably, the eontinental ice supply is mueh smaller today; Rohling et al. 

245 2008), past rapid rise is deseribed below sinee it may help inform discussions of upper bounds of 

246 future sea level rise. 

247 Average sea level rise during this more than 10,000-year period after the last ice age was 

248 9.9 to 11.9 em deeade'\ although this rise was punctuated by several shorter episodes of more 

249 rapid sea level rise. In the rapid ice-melt scenario, glaeiers and iee sheets are assumed to melt at 

250 that average rate. The meltwater term is applied as a second-order polynomial, with the average 

251 present-day iee melt rate of 1.1 cm decade'^ for 2000-2004 used as a base. This represents the 

252 sum of observed mountain glacier (Bindoff et al. 2007) and iee sheet melt (Shepherd and 

253 Wingham 2007) during this period. The rapid iee-melt scenario replaees the IPCC meltwater 

254 term with the modified meltwater term; the other three sea level terms remain unehanged. This 

^ Corrections were not made to account for reductions in glacier area over time. 
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approach does not consider how rapid ice melt might indirectly influence sea level in the New 
York region through future second-order effects including gravitational, glacial isostatic 
adjustments, and rotational terms (e.g. Mitrovica et al. 2001, 2009). 

d. Climate Projections: Extreme Events 

Based on stakeholder feedback, quantitative and qualitative projections were made using 
the extreme events definitions stakeholders currently use. For example, temperature extremes 
were defined based on specific thresholds, such as 90“F (~32°C), that the New York City 
Department of Buildings uses to define cooling requirements, whereas coastal flooding was 
defined by frequency of occurrence (Solecki et al. 2010). 

1) Quantitative projections: Coastal flood example 

The coastal flooding projections are based on changes in mean sea level, not storms. 
Projected changes in mean sea level (using the IPCC AR4-based approach) were superimposed 
onto historical data. For coastal flooding, critical thresholds for decision-making are the 1-in-lO 
year, and 1 -in- 100 year flood events (Solecki et al. 2010). The latter metric is a determinant of 
construction and environmental permitting, as well as flood insurance eligibility (Sussman and 
Major 2010). 

The 1-in-lO year event was defined using historical hourly tide data from the Battery tide 
gauge, lower Manhattan (http://tidesandcurrents.noaa.gov: for more information, see Horton and 
Rosenzweig 2010). The 1 -in- 100 year flood was analyzed using flood return period curves 
based on data provided by the U.S. Army Corps of Engineers for the Metro East Coast Regional 
Assessment (see Gomitz 2001 for details). 
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278 Because interannual variability is particularly large for rare events such as the 1-in-lO 

279 year flood, a base period of more than the standard 30 years was used. Similarly, since each year 

280 between 1962 and 1965 was drier in Central Park than the driest year between 1971 and 2000, 

281 the entire 20* century precipitation record was used for the drought analysis. More rigorous 

282 solutions for the rarest events await better predictions of interannual to multi-decadal variability, 

283 better understanding of the relationship between variability at those timescales and extreme 

284 events (see e.g. Namias 1966; Bradbury et ah, 2002) and the growing event pool of realizations 

285 with time. 

286 

287 2) Qualitative extreme event projections 

288 The question arose of how best to meet stakeholder needs when scientific understanding, 

289 data availability, and model output are incomplete; quantitative projections are unavailable for 

290 some of the important climate hazards consistently identified by infrastructure stakeholders 

291 and/or are characterized by such large uncertainties as to render quantitative projections 

292 inadvisable. Examples in the New York City region include ice storms, snowfall, lightning, 

293 intense sub-daily precipitation events, tropical storms and nor’easters. For these events, 

294 qualitative information was provided, describing only the most likely direction of change and an 

295 associated likelihood using the IPCC WGl likelihood categories (IPCC, 2007)^°. Sources of 

296 uncertainty and key historical events were also described, in order to provide stakeholders with 

Given the large impact of these extreme events on infrastructure, stakeholders requested information about 
likelihood for comparative purposes (e.g. “Which is more likely to increase in frequency? Nor’easters specifically, 
or intense precipitation events generally?”). Assignment of likelihood to generalized categories for qualitative 
extremes (based on published literature and expert judgment including peer review) was possible because 
predictions are general (e.g., direction of change),as opposed to the quantitative model-based projections. 
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297 context and the opportunity to assess sector-wide impaets of historieal extremes. 

298 

299 3. GCM Hindcasts and Observations 

300 The results of the GCM hindcasts and observational analysis deseribed in this section 

301 informed the development of the projeetion methods deseribed in seetion two. Stakeholders 

302 commonly request hindeasts and historieal analysis (see e.g. NYCDEP 2008) as they provide 

303 transparency to decision-makers who may be new to using GCM projeetions as a planning tool. 

304 

305 a. Temperature and Precipitation Trends 

306 As shown in Table 4, both the observed and modeled 20* eentury warming trends at the 

307 annual and seasonal seale are generally signifieant at the 99 pereent level. While GCM 20* 

308 eentury trends are generally approximately 50% smaller than the observed trends, it has been 

309 estimated that approximately 1/3 of New York City’s 20* eentury warming trend may be due to 

310 urban heat island effects (Gaffin et al. 2008) that are external to GCMs. Over the 1970-1999 

311 period of stronger greenhouse gas foreing, the observed annual trend was 0.2 1°C decade'\ and 

3 12 the ensemble trend was 0.1 8°C deeade'^ 

313 Modeled seasonal warming trends in the past three deeades and both annual and seasonal 

314 precipitation trends over the entire eentury for New York City generally deviate strongly from 

315 observations, eonsistent with prior results for the Northeast (see e.g. Hayhoe et al. 2007). 

316 Observed and modeled trends in temperature and precipitation at a particular location are highly 

317 dependent on internal variability, and therefore highly sensitive to the seleetion of years. For 

318 example, the 1970-1999 observed Central Park annual preeipitation trend of -1.77 em deeade'^ 

319 shifts to 0.56 em deeade'^ when the analysis is extended through 2007. This is espeeially true for 
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320 the damaging extreme events^ ^ (Christensen et al. 2007) that are often of particular interest to 

321 infrastructure managers. In coupled GCM experiments with a freely evolving climate system, 

322 anomalies associated with climate variability generally will not coincide with observations, 

323 leading to departures between observed and modeled trends (Randall et al. 2007). 

324 For stakeholders trained in analyzing recent local observations, it is challenging but 

325 important to emphasize that: 1) trends at continental and centennial timescales are often most 

326 appropriate for identifying the greenhouse gas signal and GCM performance, since 

327 (unpredictable) interannual to interdecadal variability is lower at those scales (Hegerl et al. 

328 2007); and 2) during the 2F' century, higher greenhouse gas concentrations and other radiatively 

329 important agents are expected to increase the role of the climate change signal, relative to climate 

330 variability. 

331 

332 b. Temperature and Precipitation Climatology 

333 Comparison of station data to a GCM gridbox is hindered by the spatial scale 

334 discrepancy; New York City’s low elevation, urban heat island (see e.g. Rosenzweig et al. 2006), 

335 and land sea contrasts are not captured by GCMs. As shown in Fig. 3a, the observed average 

336 annual temperature over the 1970-1999 period for New York City exceeds the GCM ensemble 

337 by 2.6°C, and is higher than all but two of the 16 GCMs. When the GCMs are contrasted with 

338 the spatially comparable NCEP Reanalysis gridbox, the annual mean temperature bias is reduced 

^ ^ Among 20* Century Central Park trends in observed extremes, only trends in eold extremes have been robust. 
For the number of days per year with minimum temperatures below freezing, both the 100-year trend of -2 days 
deeade * and the 30-year trend of -5.2 days deeade * are signifieant at the 99% level. GCM hindeasts of extreme 
events were not eondueted due to the small signal to noise ratio. 
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339 to 1.1°C. The departure of the Central Park station data from the GCM ensemble is largest in July 

340 and smallest in January, indieating that the annual temperature eyele at this location is damped in 

341 the GCMs (Fig. 3b). 

342 While Figure 3c reveals that the GCM ensemble of average annual precipitation from 

343 1970-1999 is 8% below observations for Central Park, the ensemble average lies well within the 

344 range of precipitation for New York City as a whole; GCM precipitation exceeds the LaGuardia 

345 Airport station by 9%. Most of the GCMs are able to capture the relatively even distribution of 

346 monthly precipitation throughout the year (Fig. 3d). 

347 The above analysis reveals that mean climatology departures from observations over the 

348 hindcast period are large enough to necessitate bias correction such as the delta method as part of 

349 the GCM projection approach, rather than direct use of model output. 

350 

351 c. Temperature and Precipitation Variance 

352 1) Interannual 

353 Eleven (ten) of the 16 GCMs overestimate the 1970-1999 interannual standard deviation 

354 of temperature, relative to the station data (NCEP reanalysis). The similarities between GCMs, 

355 reanalysis and station data suggest that spatial-scale discontinuities may not have a large impact 

356 on interannual temperature variance. All 16 GCMs underestimate interannual precipitation 

357 variability relative to Central Park observations, and 14 of the 16 GCMs underestimate variance 

358 relative to two other stations analyzed (Port Jervis and Bridgehampton). The large difference 

359 between the GCMs and station data suggests that spatial-scale discontinuities, likely associated 

360 with features like convective rainfall that cannot be resolved by GCMs, may be partially 

361 responsible for the relatively low modeled interannual precipitation variance. Observed 
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362 interannual temperature varianee is greatest in winter, a pattern not eaptured by seven on the 16 

363 GCMs. 

364 

365 2) High-frequency 

366 The daily distribution of observed Central Park temperature (Fig.4 a-e) and precipitation 

367 (Fig. 5) was compared to single gridbox output from 3 of the 16 GCMs used in the larger 

368 analysis. The three models were part of a subset with daily output stored in the WCRP / CMIP3 

369 repository and were selected because (of the subset) they featured the highest [Max Plank 

370 Institute for Meteorology ECHAM5/MPI-OM (MPI, Jungclaus et al. 2005) and Commonwealth 

371 Scientific and Industrial Research Organisation CSIRO-MK3.0 (CSIRO, Gordon et al. 2002), 

372 both at 1.88°lat. x 1.88“lon.] and lowest [National Aeronautics and Space Administration 

373 (NASA) / Goddard Institute for Space Studies (GISS) GISS-ER (GISS, Schmidt et al. 2006), at 

374 4°lat. X 5°lon.] resolution. Analysis was conducted on summer (June-August) daily maximum 

375 and winter (December-Eebruary) daily minimum temperature. 

376 Summer maximum temperature distribution for the region in all three GCMs is narrower 

377 than observations, and the warm tail is more poorly simulated than the cold tail. During winter, 

378 CSIRO and MPI underestimate variance relative to the station data, while the GISS GCM has 

379 excessive variance. 

380 Eigure 5 shows the number of days with precipitation exceeding 10 mm, a level of 

381 rainfall that can trigger combined sewer overflow events at vulnerable sites in New York City 

382 (PlaNYC 2008). Relative to Central Park data, all three GCMs underestimate the frequency of 

383 daily precipitation above 50 mm— a level of precipitation that can lead to widespread flooding 

384 and drainage problems including in subways (MTA 2007). 
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385 Given that precipitation in GCMs of this class and spatial resolution is highly 

386 parameterized to the gridbox spatial scale and seasonal/decadal climate timescales, departures of 

387 the distribution from observed daily station data can be expected. The low model variance at 

388 daily timescales for temperature and precipitation, and at interannual timescales for precipitation, 

389 reinforces the need for statistical downscaling approaches such as the delta method that apply 

390 monthly mean model changes to observed high frequency data. 

391 

392 d. Sea-Level Rise 

393 Sea level was also hindcast for the 20* century, based on a 1990-1999 projection relative 

394 to the 1900-1904 base period. The ensemble average hindcast is a rise of 18 cm, while the 

395 observed increase at the Battery is 25 cm. The five-year average local elevation term in the 

396 models meanders through time, frequently with an amplitude of 2-3 cm, with a maximum range 

397 over the century of approximately 7 cm, suggesting decadal variability (primarily in the local 

398 elevation term) and spatial resolution may explain the discrepancy between models and 

399 observations. 

400 

401 4. Future Projections 

402 a. Mean Temperature and Precipitation 
1 2 

In this calculation, the land subsidence term was identical to that used for the 2fi* century projections. The same 
surface mass balance coefficients used by the IPCC, based on global average temperature changes over a 1961-2003 
baseline were used for the 1900-1904 base period, which likely leads to a slight overestimate of the meltwater here. 
The effect is negligible though as the meltwater term is a minor contributor to the overall 20* century sea level rise. 
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403 1) Annual 

404 Table 5 shows the projected changes in temperature and precipitation for the 30-year 

405 periods centered around the 2020s, 2050s, and 2080s relative to the baseline period. The values 

406 shown are the central range (middle 67%) of the projected model-based changes. 

407 Figure 6 expands upon the information presented in Table 5 in three ways. First, 

408 inclusion of observed data since 1900 provides context on how the scale of projected changes 

409 associated with forcing from greenhouse gases and other radiatively important agents compares 

410 to historical variations and trends. Secondly, tabulating high and low projections across all 48 

411 simulations provides a broader range of possible outcomes, which some stakeholders requested 

412 (New York City Climate Change Adaptation Task Force, 2008-2009). Finally, ensemble 

413 averaging of results by emissions scenario as they evolve over time is informative to 

414 stakeholders involved in greenhouse gas mitigation (and adaptation), since it reveals the large 

415 system inertia: not until the 2030s and 2040s do the B1 scenario projections begin to diverge 

416 from A2 and AIB, but thereafter it diverges rapidly. Thus, a delay in greenhouse gas mitigation 

417 activities greatly increases the risk of severe long-term climate change consequences, despite 

418 apparent similarity in the near-term outlook. 

419 While the precise numbers in Table 5 and Figure 6 should not be emphasized due to high 

420 uncertainty and the smoothing effects of ensemble averaging, the stakeholder sees that in the 

421 New York Metropolitan Region: 1) mean temperatures and sea levels are projected to increase in 

422 all simulations this century, at rates exceeding those experienced in the 20* century; 2) while 

423 precipitation is projected to increase slightly in most simulations, the multi-year precipitation 

424 range experienced in the past century due to climate variability exceeds the 2F' century climate 
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425 change signal^^; and 3) climate projection uncertainties grow throughout the 21®' century, in step 

426 with uncertainties regarding future emissions and the climate system response. 

427 

428 2) Seasonal 

429 Warming in the New York City region is of similar magnitude for all seasons in the 

430 GCMs, although seasonal projections are characterized by larger uncertainties than annual 

431 projections (Fig. 7a). Since interannual temperature variability is smallest in summer this 

432 suggests the summer warming may produce the largest departures from historical experience. 

433 Some impacts and vulnerabilities are also amplified by high temperatures. Energy demand in 

434 New York City is highly sensitive to temperature during heat waves, due especially to increased 

435 reliance on air conditioning. This increased demand can lead to elevated risk of power shortages 

436 and failures at a time when vulnerable populations are exposed to high heat stress and air 

437 pollution (Kinney et al. 2001; Kalkstein 1995; Hill and Goldberg 2001; Hogrefe et al. 2004). 

438 GCMs tend to distribute much of the additional precipitation during the winter months 

439 (Fig. 7b), when water supply tends to be relatively high and demand relatively low (NYCDEP 

440 2008). During September and October, a time of relatively high drought risk, total precipitation 

441 is projected to decrease slightly in many models. 

442 

The projection lines in Figure 6 depict the ‘predictable’ anthropogenic forcing component, while capturing some 
of the uncertainty associated with greenhouse gas concentrations and climate sensitivity at specific points in time. 
Because decadal variability is unpredictable in the Northeast, it was not included in the time-specific projection 
portion of the figure. It was however emphasized to stakeholders that while interannual variability appears greatly 
reduced in the projection portion of the figure, the observed portion (black line) reflects the type of unpredictable 
variations that have been experienced in the past and will likely exist on top of the mean change signal in the future. 
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443 b. Sea Level Rise 

444 Addition of the two regional eomponents leads to higher sea level rise projeetions for the 

445 region than the global average (by '-'15 em for end-of-century projections; Meehl et al. 2007b; 

446 Peltier 2001). This is due both to land subsidence and higher sea level rise along the northeast 

447 U.S. coast, the latter largely due to geostrophic constraints associated with projected weakening 

448 of the Gulf Stream (Yin et al. 2009) in many GCMs (Meehl et al. 2007b). 

449 As shown in Table 6, the rapid ice melt scenario projections diverge from the IPCC- 

450 based approach as the century progresses. The 2100 value of up to ~2 meters associated with 

451 this scenario (not shown) is generally consistent with other recent results that roughly constrain 

452 sea level rise globally (see e.g., Pfeffer et al. 2008; Rahmstorf, 2007; Horton et al. 2008; Grinsted 

453 2009; Rignot and Cazenave 2009) and regionally (see Yin et al. 2009; Hu et al. 2009) to between 

454 --Im and -'2m. The consistency with other studies supports the usefulness of -'2m as a high end 

455 for a risk-averse approach to century-scale infrastructure investments including bridges and 

456 tunnels, rail lines, and water infrastructure. 

457 At the request of agencies that manage some of these long-term investments, two 

458 presentations were given to technical staff specifically describing the rapid ice melt methodology 

459 and projections. While these and other stakeholders wanted to know the probability of the rapid 

460 ice melt scenario relative to the IPCC-based method, it was emphasized that such probability 

461 statements are not possible given current scientific understanding. 

462 

463 c. Extreme Events 

464 1) STAKEHOLDER PROJECTIONS BASED ON THE DELTA METHOD 

465 Table 7 shows projected changes in the frequency of heat waves, cold events, intense 
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466 precipitation, and coastal flooding in the New York City region. The baseline average number of 

467 extreme events per year is shown, along with the central range (middle 67%) of the projections. 

468 Because the distribution of extreme events around the (shifting) mean could also change while 

469 mean temperature, precipitation, and sea level rise shift, stakeholders were strongly encouraged 

470 to focus only on the direction and relative magnitudes of the extreme event changes in Table 7. 

471 The key finding for most stakeholders is the extent to which mean shifts alone can 

472 produce dramatic changes in the frequency of extreme events, such as heat events and coastal 

473 storm surges. Based on the central range, the number of days per year over 90 “F is projected to 

474 increase by a factor of approximately three by the 2080s. The IPCC-based sea level rise 

475 projections alone, without any changes in the historical storm climatology and surge levels, lead 

476 to a more than threefold increase in the frequency of the baseline 1-in-lO year coastal flood event 

477 by the 2080s. 

478 In contrast to relatively homogeneous mean climate changes, it was emphasized to 

479 stakeholders that absolute extreme event projections like days below freezing and days with 

480 more than one inch of precipitation vary dramatically throughout the metropolitan region, since 

481 they depend for example on microclimates associated with the urban heat island and proximity to 

482 the coast. Similarly, maps were generated for stakeholders to show that the surge heights for the 

483 open estuary at the Battery are higher than corresponding heights in more protected riverine 

484 settings. 

485 It was emphasized to stakeholders that due to large interannual variability in extremes, 

486 even as the climate change signal strengthens, years with relatively few extreme heat events 

487 (relative to today’s climatology) will occur. For example. Central Park’s temperatures in 2004 

488 only exceeded 90°F (~32°C) twice. The delta method suggests that not until the middle of this 
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489 century would such a relatively cool summer (as 2004) feature more days above 90°F (~32°C) 

490 than are typically experienced today. 

491 High year-to-year extreme event variability may already give some stakeholders a 

492 framework for assessing sector-specific climate change impacts; even if climate adaptation 

493 strategies for extremes are not already in place, short-term benefits may be evident to planners. 

494 For example, Central Park in 2010 experienced temperatures of higher than 90°F (-'32°C), on 32 

495 different days, which is consistent with projections for a typical year around mid-century. This 

496 suggests that some of the infrastructure impacts of extreme heat (such as voltage fluctuations 

497 along sagging power lines and increased strain on transportation materials including rails and 

498 asphalt; Horton and Rosenzweig 2010) may have been experienced in 2010 to an extent that may 

499 become typical by mid-century. However, adaptation strategies designed for an extreme year 

500 today (such as a fixed level of mandatory energy use reductions and a fixed level of reductions of 

501 train speeds) may be inadequate or unpalatable in the future due to the increase in frequency, 

502 duration and intensity of extreme heat (for example) associated with climate change (see e.g., 

503 Meehl et al. 2009; Tebaldi et al. 2006; Meehl and Tebaldi 2004). 

504 

505 2) GCM changes in intra-annual distributions 

506 Since high frequency events are not well-simulated in GCMs, the results described here 

507 were not included in the New York City adaptation assessment; they are explored here as an 

508 exercise, since there is the possibility of distributional changes in the future. The daily 

509 distribution of: a) maximum temperatures^^ in summer (JJA), and b) minimum temperatures in 

Precipitation was excluded, based on the preliminary analysis of hindcast daily precipitation described in section 
3d. 
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winter (DJF) are analyzed in the three GCMs described earlier (CSIRO, GISS and ECHAM5; 
section 3d), both for the 1980-1999 hindcast and the 2080-2099 AIB experiment. 

The results indicate that GCM temperature changes in the region in some cases do reflect 
more than a shifting mean. The intra-aimual standard deviation of winter minima decreases in 
all three GCMs (in two cases by approximately 10 %), while summer standard deviation changes 
are negligible. One tail of a season’s distribution can be more affected than the other; as shown 
in Fig. 8 for CSIRO, the winter minimum changes are more pronounced on anomalously cold 
days than anomalously warm days. All 3 GCMs show a larger shift in the coldest 1% of the 
distribution than the highest 1%. This asymmetry at the 1% tails is most pronounced in CSIRO, 
where the future coldest 1% event occurs 8 times more often in the baseline, while the baseline 
warmest 1% event occurs three times more often in the future. 

d. Comparison of GCM gridbox based-projections to other downscaling methods 

The GCM grid box results used for the New York assessment were compared to 
statistically downscaled results from Bias-Corrected and Spatially -Downscaled (BCSD) Climate 
Projections at 1/8 degree resolution derived from the World Climate Research Programme's 
(WCRP's) Coupled Model Intercomparison Project Phase 3 (CMIP3) multi-model dataset. The 
BCSD projections are available at: http://gdo-dcp.ucllnl.org/downscaled cmip3 projections/ 
(Maurer, 2007). Results were also compared to simulations from four parings of GCMs and 
RCMs (Table 8) contributing to the North American Regional Climate Change Assessment 
Program (NARCCAP; Meams et al. 2009). Comparison of the three methods is limited to the 

As calculated separately for each year and then averaged across the 20 years, to minimize the role of interannual 
variability. 
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531 2050s timeslice under the A2 emissions seenario relative to the 1970-1999 baseline, sinee 

532 NARCCAP projections are not available for other emissions scenarios or time periods. The 

533 comparison focuses on projections rather than validation, since the BCSD methodology by 

534 definition includes bias correction whereby the baseline GCM outputs are adjusted to match the 

535 observed mean and variance. Preliminary analysis of NARCCAP results indicates that these 

536 simulations, like GCM projections, require bias correction. 

537 The ensemble mean changes for the GCM gridbox, BCSD, and RCM approaches differ 

538 from each other by no more than .3“C for temperature and 3% for precipitation. The intermodel 

539 temperature range is slightly larger for the GCM gridbox approach than BCSD, while the 

540 opposite is the case for precipitation. The four RCM simulations perhaps not surprisingly feature 

541 a smaller intermodel range than the 16 ensemble members for the GCM gridbox and BCSD 

542 approaches. 

543 The number of days above 90 °F was evaluated as a measure of extremes events. The 

544 delta method applied to the GCM gridbox and BCSD^^ produce virtually identical results 

545 (increases of approximately 185 and 180 percent respectively in the number of days above 90°F). 

546 When actual daily values from RCMs are used the increase is approximately 170 percent. When 

547 the delta method from the RCMs is applied to the observations, the increase is approximately 

548 195 percent. 

549 For mean changes and the daily extreme metric assessed here, BCSD and the four RCMs 

550 offer comparable results to the single gridbox GCM approach in the New York Metropolitan 

551 Region. Future research will assess how statistical and dynamic downscaling perform in more 

552 specialized contexts tailored to unique stakeholder needs that are beyond the scope of New York 

At the time of analysis, BCSD is only available at monthly resolution. 
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553 City initial assessment. For example, reservoir managers eoneerned with water turbidity might 

554 desire information about sequenees of days with intense preeipitation during partieular times of 

555 year. Future researeh will also explore the pros and eons of projeetions that ineorporate highly 

556 uneertain modeled ehanges in interannual variance through time^’. 

557 

558 5. Conclusions and Recommendations for Future Work 

559 A framework for climate hazard assessment geared towards adaptation planning and 

560 decision support is described. This single GCM gridbox, delta method-based approach, designed 

561 for cities and regions smaller than typical GCM gridbox sizes that face resource and time 

562 constraints, achieves comparable results in the New York Metropolitan Region to other 

563 statistically and dynamically downscaled products. When applied to high frequency historical 

564 data, long-term mean monthly climate changes (which GCMs are expected to simulate more 

565 realistically for point locations than other features such as actual long-term mean climate or high 

566 frequency statistics) yield dramatic changes in the frequency of stakeholder-relevant climate 

567 hazards such as coastal flooding and heat events. While the precise projections should not be 

568 emphasized given the uncertainties, they are of sufficient magnitude relative to the historical 

569 hazard profile to justify development and initial prioritization of adaptation strategies. This 

570 process is now well underway in the New Y ork Metropolitan Region. 

571 When climate model results for the New York Metropolitan Region are used only for the 

572 calculation of monthly climate change factors based on the differences and ratios between 30- 

Preliminary analysis reveals that over the New York Metropolitan Region gridbox a slight majority of the GCMs 
show increasing interannual variance of monthly T and P, while a large majority of the BCSD and NARCCAP RCM 
projections do. 
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573 year future timesliees and a 30-year baseline period, three generalized findings follow. First, 

574 using multiple ensemble members from the same GCM provides little additional information, 

575 sinee the 30-year average intramodel ranges are smaller than the eomparable inter-model range. 

576 Seeond, the spatial pattern of climate change factors in many regions (including New York City) 

577 is sufficiently homogeneous — relative to the intermodel range — to justify use of climate 

578 change factors from a single overlying GCM gridbox. Finally, for these metrics, newer 

579 statistically (BCSD) and dynamically (four NARCCAP RCMs) downscaled products provide 

580 comparable results to the GCM single gridbox output used by the NPCC. 

581 The checklist in Table 2 provides a series of questions to help inform the selection of the 

582 most appropriate climate hazard assessment and projection methods. For example, the delta 

583 method is more justified when: 1) robust, long-term historical statistics are available, and 2) 

584 evidence of how modes of interannual and interdecadal variability and their local teleconnections 

585 will change with climate change is inconclusive. Both these criteria are met in the New York 

586 City Metropolitan Region. In contrast, more complex applications (than the delta method) of 

587 statistically and dynamically downscaled products especially may be more appropriate when 

588 spatially continuous projections are needed over larger regions with complex topography. For 

589 example, where a large mountain range is associated with a strong precipitation gradient at sub- 

590 GCM gridbox scales, percentage changes in precipitation might also be expected to be more 

591 spatially heterogeneous than in the New York Metropolitan Region. 

592 Extreme event projections, so frequently sought by stakeholders for impact analysis, will 

593 likely improve as statistical and dynamical downscaling evolve. RCMs especially hold promise 

594 for assessing how ‘slow’ variations associated with climate change and variability will affect the 

595 future distribution of ‘fast’ extremes like subdaily rainfall events. Nevertheless, translating RCM 
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596 simulations into stakeholder-relevant projeetions requires many of the same adjustments and 

597 caveats described here for GCMs (such as bias correction). Statistical downscaling techniques 

598 also hold promise as well for the simulation of extremes (non-stationarity not withstanding), to 

599 the extent that predictor variables are well simulated by GCMs and linkable to policy relevant 

600 local climate variables. Projections of extremes will also benefit from improved estimates of 

601 historical extremes (such as the 1-in-lOO year drought and coastal flood) as long-term tree ring 

602 and sediment records (for example) are increasingly utilized. 

603 There is also a need for improved simulation of climate variability at interannual to 

604 decadal scales, as this is the time horizon for investment decisions and infrastructure lifetime in 

605 many sectors, including telecommunications (NPCC, 2010). The limits to such predictability are 

606 beginning to be explored in Coupled Model Intercomparison Project (CMIP5) experiments 

607 initialized with observed ocean data, but this is a long-term research issue. 

608 An absence of local climate projections need not preclude consideration of adaptation. 

609 For many locales, climate changes in other regions may rival the importance of local changes by 

610 influencing migration, trade, and ecosystem and human health, for example. Furthermore, some 

611 hazards such as drought are often regional phenomena, with multi-state policy implications (such 

612 as water-sharing agreements). Finally, since climate vulnerability depends on many non-climatic 

613 factors (such as poverty), some adaptation strategies (such as poverty-reduction measures) can be 

614 commenced in advance of climate projections. 

615 Monitoring of climate indicators should be encouraged since it reduces uncertainties and 

616 leads to refined projections. Locally, sustained high temporal resolution observation networks 

617 can provide needed microclimatic information, including spatial and temporal variation in 

618 extreme events such as convective rainfall and storm surge propagation. At the global scale. 
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619 monitoring of polar ice sheets and global sea level will improve understanding of sea level rise. 

620 Periodic assessments of evolving elimate, impaets and adaptation seienee will support 

621 flexible/recursive adaptation strategies that minimize the impact of climate hazards while 

622 maximizing soeietal benefits. 

623 
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915 FIGURE CAPTIONS 

916 Figure 1; Satellite map of the New York Metropolitan Region. Shown on the map are the 

917 Central Park weather station (cirele) and The Battery tide gauge (triangle). Source: ESRl World 

918 Imagery. 

919 

920 Figure 2: a) Temperature change (°C) and b) precipitation change (%) for the 2080s timeslice 

921 relative to the 1970-1999 model baseline, AIB emissions scenario and 16 GCM ensemble mean. 

922 

923 Figure 3: a) Mean annual temperature for the New York City region, (°C), 1970-1999 in each of 

924 the 16 GCMs, GCM ensemble. Central Park station data and Reanalysis (see methods section for 

925 more information). Also shown as hash marks is the interannual standard deviation about the 

926 mean for each of the 19 products, b) monthly mean temperature for the New York City region, 

927 (°C), 1970-1999. The two observed products, the GCM ensemble average, and four points in the 

928 GCM distribution (lowest, 17th percentile, 83rd percentile, and highest) are shown, c) Mean 

929 annual precipitation for the New York City region, (cm), 1970-1999 in each of the 16 GCMs, 

930 GCM ensemble, and Central Park observations. Also shown as hash marks is the interannual 

931 standard deviation about the mean for each of the 18 products, d) monthly mean precipitation for 

932 the New York City region, (cm), 1970-1999. Central park observations, the GCM ensemble 

933 average, and four points in the GCM distribution (lowest, 17th percentile, 83rd percentile, and 

934 highest) are shown. 

935 

936 Figure 4: Daily distribution (number of days per year) of: a) all-year mean, b) summer (June- 

937 August) maximum, and c) winter (December-February) minimum temperature anomalies (°C), 
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938 1980-1999 for Central Park observations (blaek line) and three GCMs (CSIRO, GISS, and MPl 

939 ECHAM5). 

940 

941 Figure 5: Daily distribution (number of days per year) of preeipitation (mm), 1980-1999 for 

942 Central Park observations (blaek line) and three GCMs (CSIRO, GISS, and MPl ECHAM5). The 

943 first bin, eontaining less than 10 mm, is not shown. 

944 

945 Figure 6: Combined observed (blaek line) and projeeted: a) temperature (°C) and b) annual 

946 preeipitation (mm). Projeeted model ehanges through time are applied to the observed historieal 

947 data. The three thiek lines (red, green, and blue) show the ensemble average for eaeh emissions 

948 seenario aeross the 16 GCMs. Shading shows the eentral 67 % range aeross the 16 GCMs and 3 

949 emissions seenarios. The bottom and top lines, respeetively, show eaeh year’s minimum and 

950 maximum projeetions aeross the suite of simulations. A ten -year filter has been applied to the 

951 observed data and model output. The dotted area between 2003 and 2015 represents the period 

952 that is not eovered due to the smoothing proeedure. 

953 

954 Figure 7: Seasonal a) temperature ehange (°C) and b) preeipitation ehange (%) projeetions, 

955 relative to the 1970-1999 model baseline, based on 16 GCMs and 3 emissions seenarios. The 

956 maximum and minimum are shown as blaek horizontal lines; the eentral 67% of values are 

957 boxed, and the median is the thiek line inside the boxes. 

958 

959 Figure 8; Daily distribution (number of days per year) of winter (Deeember-February) minimum 

960 temperature anomalies (°C), for the New York Metropolitan Region in the CSIRO GCM. Blaek 
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961 line, 1980-1999 hindcast; dotted line, 2080-2099 AIB seenario. 
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962 TABLE 1. Acronym, host center, atmosphere and oeean grid box resolution, and referenee for the 16 GCMs used in the analysis. 


Model 

Institution 

Atmospherie 

Oeeanie 

Referenees 

Aeronym 


Resolution 

Resolution 




(Eat X Eon ) 

(Eat X Eon) 


BCCR 

Bjerknes Center for Climate Researeh, Norway 

1.9 X 1.9 

0.5 to 1.5 X 1.5 

Furevik et ah, 2003 

CCSM 

National Center for Atmospherie Researeh, USA 

1.4 X 1.4 

0.3 to 1.0 X 1.0 

Collins et ah, 2006 

CGCM 

Canadian Center for Climate Modeling and Analysis, 
Canada 

2.8 X 2.8 

1.9 X 1.9 

Flato 2005 

CNRM 

National Weather Researeh Center, METEO-FRANCE, 
Franee 

2.8 X 2.8 

0.5 to 2.0 X 2.0 

Terray et ah, 1998 

CSIRO 

CSIRO Atmospherie Researeh, Australia 

1.9 X 1.9 

0.8 X 1.9 

Gordon et ah, 2002 

ECHAM5 

Max Planek Institute for Meteorology, Germany 

1.9 X 1.9 

1.5 X 1.5 

Jungelaus et ah, 2005 

ECHO-G 

Meteorologieal Institute of the University of Bonn, 
Germany 

3.75x3.75 

0.5 to 2.8 X 2.8 

Min et ah, 2005 

GEDE-CM2.0 

Geophysieal Fluid Dynamies Eaboratory, USA 

2.0 X 2.5 

0.3 to 1.0 X 1.0 

Delworth et ah, 2006 

GEDE-CM2.1 

Geophysieal Fluid Dynamies Eaboratory, USA 

2.0 X 2.5 

0.3 to 1.0 X 1.0 

Delworth et ah, 2006 
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GISS 

NASA Goddard Institute for Space Studies 

4.0 X 5.0 

4.0 X 5.0 

Schmidt et ah, 2006 





Volodin and Diansky, 

INMCM 

Institute for Numerical Mathematics, Russia 

4.0 X 5.0 

2.0 X 2.5 

2004 

IPSL 

Pierre Simon Laplace Institute, France 

2.5x3.75 

2.0 X 2.0 

Marti, 2005 

MIROC 

Frontier Research Center for Global Change, Japan 

2.8 X 2.8 

0.5 to 1.4 X 1.4 

K-1 Developers, 2004 


Meteorological Research Institute, Japan 



Yuikimoto and No da. 

MRI 


2.8 X 2.8 

0.5 to 2.0 X 2.5 

2003 

PCM 

National Center for Atmospheric Research, USA 

2.8 X 2.8 

0.5 to 0.7 X 1.1 

Washington et ah, 2000 

UKMO- 

Hadley Center for Climate Prediction, Met Office, UK 

2.5x3.75 

1.25 X 1.25 

Johns et ah, 2006 


HadCM3 
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964 Table 2. Checklist of questions to inform selection of climate hazard assessment and projection methods 

Question Possible implication for choice of method, plus NYC context 

1. Are high quality historical data available for a long time period? When little high-quality historical climate data are available, 

options for projections are extremely limited. Records of at least 
several decades are needed to sample the range of natural 
variability. As RCMs continue to improve, use of raw outputs 
from RCMs may increasingly be used in such regions, since bias 
correction and statistical approaches are not feasible without 
historical climate data. This was not an issue in data-rich NYC. 

2. Are projections needed for the entire century? If yes, this may preclude RCMs due to computational expense. 

This was an important consideration for NYC, since some sectors 
such as telecommunications were focused on the 2020s timeslice, 
while others such as Port Authority of New York and New Jersey 
manage infrastructure expected to last until 2100. 

2. Are multiple emissions scenarios needed, for example to If yes, RCMs may not be the best approach, since computational 

emphasize how mitigation can compliment adaptation? expense generally precludes the use of more than 1-2 scenarios. 
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4. Are a large group of GCMs and initializations required, in order 
to sample a broad range of global climate sensitivities and estimates 
of within-GCM variability, respectively? 

5. What climate variables are needed, and are they available at the 
necessary spatial and temporal resolution within public climate 


model archives? 


This was an important consideration in New York City, sinee the 
adaptation effort was part of a broader sustainability effort 
(PlaNYC) that embraeed greenhouse gas mitigation. 

If yes, RCMs may not be the best approach, since eomputational 
expense generally preeludes the use of more than a few GCMs or 
GCM initializations per RCM. New York City stakeholders 
expressed interest in the full range of GCM sensitivities. 

In NYC, relatively few variables were needed and subdaily 
information was not required. Additional variable needs at 
subdaily resolution might argue for the use of RCM archives sueh 
as NARCCAP as they eontinue to be populated, instead of 
arehives such as the first generation of BCSD (monthly 
temperature and preeipitation only). While use of publie climate 
model arehives minimizes cost and time, even archived outputs 
generally require at least some bias and/or seale correetion and 
post-proeessing for stakeholder applieability. 
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6. What level of resources are available, and in what time frame is 
the information needed? 


7. Are projections needed for a single in-depth sectoral application 
and variable in one municipality, or does a large multisectoral and 
pan-regional group of stakeholders need a coordinated set of 
scenarios covering a series of standard variables? 


8. Are high-frequency climate inputs that are continuous in time and 


Region and question speeifie tailored downsealing efforts, as 
opposed to use of arehived downsealed produets, may not be 
possible when resourees and time are limited. While NY C had 
substantial resourees available, the short time frame ('-8 months) 
preeluded developing new tailored downsealing. 

In tailored statistieal downsealing the method is optimized to the 
partieular loeation and/or variable. When many variables and a 
larger region are ineluded, no single optimization method will 
generally be best for all variables and loeations, potentially 
leading to ineonsisteneies in either methods or projeetions aeross 
variables and loeations. In NYC, the initial emphasis was on 
generating a eommon denominator of eonsistent seenarios based 
on eonsistent methods (the delta method) to faeilitate eoordination 
aeross 40 stakeholder entities. 

If an impaets model is to be run with elimate outputs, the range of 
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space required, such as for input into an impacts model? (e.g., a 
hydrological model to assess turbidity) 


9. Is the region ’s climate characterized by large spatial 
heterogeneity? 

10. Are modes of variability important and predictable? 


climate and impact results (rather than just the ‘delta’ mean) will 
likely be of interest, which may argue for a downscaling 
technique that allows variance to change, such as BCSD. 
Statistical downscaling techniques that include weather generators 
(such as SDSM) may be desirable to create a long record at the 
needed resolution that includes a range of extreme outcomes for 
planning purposes. The larger the continuous geographic domain 
(e.g., a large watershed) the greater the need for caution regarding 
weather generator treatment of spatio-temporal correlation. 

While impacts modeling was not the initial thrust of the NYC 
CCATF effort, climate scenarios for impact modeling are being 
developed for specific sectors (e.g., NYCDEP, 2008). 

If not, applying the delta method to a single GCM gridbox may be 
justifiable for many applications, as it was in NYC. 

If not, the use of 30-year time slices (and the delta method) that 
emphasize the signal of greenhouse gases and other radiatively 
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important agents should be emphasized, as it was in NYC. 


965 

966 
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967 TABLE 3. NCAR CCSM climatology of available simulations and CCSM ensemble for the 

968 gridbox covering New York City: a) 1970-1999 hindeast; b) AIB 2080s (2070 to 2099 average) 

969 relative to the same-simulation 1970 to 1999 hindeast. 

970 a) 



1970 - 1999 

1970 - 1999 


Mean 

Mean 


temperature 

preeipitation 


(°C) 

(cm) 

CCSM Runl 

9.38 

98.03 

CCSM Run 2 

9.27 

91.88 

CCSM Run 3 

9.67 

92.08 

CCSM Run 5 

9.42 

94.87 

CCSM Run 6 

9.64 

95.22 

CCSM Run 7 

9.64 

91.30 

CCSM Run 9 

9.68 

94.69 

CCSM 

Ensemble 

9.53 

94.10 

b) 


2080s AIB 

2080s AIB 


Temperature 

Preeipitation 


change (°C) 

change (%) 

CCSM Runl 

3.44 

2.81 
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CCSM Run 2 

3.32 

10.15 

CCSM Run 3 

3.03 

12.44 

CCSM Run 5 

3.24 

9.75 

CCSM Run 6 

2.75 

9.56 

CCSM Run 7 

2.96 

12.03 

CCSM Run 9 

3.01 

10.36 

CCSM 

Ensemble 

3.11 

9.52 
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973 TABLE 4. Annual and seasonal temperature (a,b,( °C deeade'^)) and preeipitation (c,d, (em 

974 deeade'^) ) trends, and 20th eentury (a,e) and 1970-1999 (b,d). Shown are observed Central Park 

975 station data, the 16 GCM ensemble, and four points on the GCM distribution (lowest, 17th 

976 pereentile, 83rd pereentile, and highest). 

977 a) 


20th 

century* 

Min 

16% 

83% 

Max 

Ensemble 

Observed 

Annual 

-0.03 

0.02 

0.12 

0.17 

0.07" 

0.15" 

DJF 

-0.04 

0.02 

0.16 

0.19 

0.08 

0.20 

MAM 

-0.05 

-0.02 

0.12 

0.25 

0.06 

0.18 

JJA 

-0.02 

0.03 

0.11 

0.15 

0.07 

0.12 

SON 

0.00 

0.03 

0.15 

0.18 

0.09 

0.08 

b) 

1970 - 1999 

Min 

16% 

83% 

Max 

Ensemble 

Observed 

Annual 

-0.11 

0.10 

0.28 

0.39 

0.18" 

0.21 

DJF 

-0.47 

-0.05 

0.35 

0.51 

0.11 

0.76 

MAM 

-0.36 

-0.15 

0.41 

0.74 

0.14 

0.10 

JJA 

-0.01 

0.13 

0.29 

0.44 

0.20 

0.05 

SON 

-0.06 

0.13 

0.50 

0.70 

0.29 

-0.03 


980 

981 

982 
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983 


c) 


20th 

eentury* 

Min 

16% 

83% 

Max 

Ensemble Observed 

Annual 

-1.22 

-0.22 

0.66 

0.76 

0.16 

1.60 

DJF 

-0.23 

-0.18 

0.27 

0.78 

0.05 

0.27 

MAM 

-0.27 

-0.13 

0.28 

0.39 

0.10 

0.90 

JJA 

-0.69 

-0.39 

0.22 

0.35 

-0.07 

-0.09 

SON 

-0.25 

-0.08 

0.32 

0.46 

0.10 

0.61 

d) 

1970 - 1999 

Min 

16% 

83% 

Max 

Ensemble Observed 

Annual 

-3.52 

0.02 

2.05 

5.73 

0.87 

-1.77 

DJF 

-3.21 

-0.19 

1.48 

2.94 

0.48 

-0.48 

MAM 

-2.33 

-1.37 

1.05 

1.98 

-0.08 

1.55 

JJA 

-2.08 

-1.33 

1.19 

1.75 

-0.03 

-1.51 

SON 

-1.72 

-0.55 

1.89 

2.93 

0.48 

-1.72 


986 

987 

988 * Only 15 GCMs were available for the 1900-1999 hindeast. 

989 ** Trend is signilieant at the 99% level. 

990 

991 

992 
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993 TABLE 5. Mean annual changes in temperature and precipitation for New York City.* 



2020s 

2050s 

2080s 

Air temperature ** 

+ 0.8 to 1.7° C 

+ 1.7 to 2.8° C 

+ 2.2 to 4.2° C 

Precipitation ** 

+ 0 to 5 % 

+ 0to 10% 

+ 5 to 10% 


994 

995 * Based on 16 GCMs and 3 emissions scenarios. 

996 ** Shown is the central range (middle 67%) of values from model-based distributions; 

997 temperatures ranges are rounded to the nearest tenth of a degree and precipitation to the nearest 

998 5%. 

999 
1000 
1001 
1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 
1011 
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1012 TABLE 6. Sea level rise projections for New York City^ 



2020s 

2050s 

2080s 

IPCC-based"’ 

+ 5 to 13 cm 

+ 18 to 30 cm 

+ 30 to 54 cm 

Rapid ice-melt scenario'^’‘^ 

13 to 25 cm 

48 to 74 cm 

104 to 140 cm 


1013 

1014 ^Based on 7 GCMs and 3 emissions scenarios. 

1015 ^ Shown is the central range (middle 67%) of values from model-based distributions rounded to 

1016 the nearest cm. 

1017 Rapid ice-melt scenario is based on recent rates of ice melt in the Greenland and West Antarctic 

1018 lee sheets and paleoclimate studies. See text for details. 

1019 

1020 
1021 
1022 

1023 

1024 

1025 

1026 

1027 

1028 

1029 

1030 

1031 


59 



1032 TABLE 7. Extreme Events Projeetions. 



Extreme Event 

Baseline 

2020s 

2050s 

2080s 



(1971- 2000) 





# of days per year with 





Heat & cold 

maximum temperature 





events ^ 

exceeding: 






90°E (~32°C) 
100°E (~38°C) 

14 

0.4^ 

23 to 29 
0.6 to 1 

29 to 45 
1 to 4 

37 to 64 
2 to 9 


# of heat waves per year® 

2 

3 to 4 

4 to 6 

5 to 8 


Average duration (in days) 

4 

4 to 5 

5 

5 to 7 


# of days per year with 






minimum temperature at 
or below 32°E (0°C) 

72 

53 to 61 

45 to 54 

36 to 49 

Coastal 
floods & 
storms'^’® 

1-in-lO yr flood to reoccur, 
on average 
Elood heights (in m) 
associated with 1-in-lO yr 
flood 

'-once every 
10 yrs 

1.9 

-once every 
8 to 10 yrs 

2.0 to 2.1 

-once every 
3 to 6 yrs 

2.1 to 2.2 

-once every 
1 to 3 yrs 

2.3 to 2.5 


1 -in- 100 yr flood to 

'-once every 

-once every 

-once every 

-once every 


reoccur, on average 

100 yrs 

65 to 80 yrs 35 to 55 yrs 

15 to 35 yrs 
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Flood heights (in m) 

associated with 1 -in- 100 yr 2.6 2.7 to 2.7 2.8 to 2.9 2.9 to 3.2 

flood 

1033 ‘‘ Shown is the central range (middle 67%) of values from model-based distributions based on 16 GCMs 

1034 and 3 emissions scenarios. 

1035 ’’ Decimal places shown for values less than 1 (and for all flood heights) 

1036 Defined as 3 or more consecutive days with maximum temperature exceeding 90°F (~32°C). 

1037 Does not include the rapid ice-melt scenario. 

1038 ® Shown is the central range (middle 67%) of values from model-based distributions based on 7 

1039 GCMs and 3 emissions scenarios. 
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1040 TABLE 8. Global climate model and regional elimate model pairings used from NARCCAP 


Global Climate Model 
Driver 

Regional Climate Model 

Combination 

RCM 

Referenee 

Geophysieal Fluid 
Dynamies Laboratory 
(GFDL) 

Regional Climate Model Version 3 
(RCM3) 

RCM3 + 
GFDL 

Pal et ah, 
2007 

Third Generation 
Coupled Climate Model 
(CGCM3) 

Regional Climate Model Version 3 
(RCM3) 

RCM3 + 
CGCM3 

Pal et ah, 
2007 

Third Generation 
Coupled Climate Model 
(CGCM3) 

Canadian Regional Climate Model 
(CRCM) 

CRCM + 
CGCM3 

Caya and 
Laprise, 1999 

Hadley Centre Coupled 
Model, Version 3 
(HadCM3) 

Hadley Regional Model 3 / 
Providing Regional Climates for 
Impaets Studies (HRM3) 

HRM3 + 
HadCM3 

Jones et ah, 
2004 


1041 
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1042 


1043 



Bronx 


Manhattan 


Queens 


Brooklyn 


Staten Island 


Legend 

• Central Park 
▲ The Battery 


1044 


1045 

1046 

1047 


Figure 1; Locations of the Central Park weather station (circle), and The Battery tide gauge 
(triangle), and the 5 boroughs of New York City. Source; ESRI World Imagery. 
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1049 

1050 

1051 

1052 

1053 

1054 

1055 

1056 

1057 

1058 

1059 

1060 

1061 

1062 

1063 

1064 

1065 


a) 




2 3 4 -2 2 6 10 

Annual temperature change (°C) Annual precipitation change (%) 


Figure 2: a) Temperature ehange (°C) and b) preeipitation ehange (%) for the 2080s timesliee 
relative to the 1970-1999 model baseline, AIB emissions seenario and 16 GCM ensemble mean. 
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1066 

1067 

1068 

1069 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1078 

1079 

1080 

1081 

1082 

1083 






Figure 3; a) Mean annual temperature for the New York City region, (“C), 1970-1999 in eaeh of 
the 16 GCMs, GCM ensemble, Central Park station data and Reanalysis (see methods section for 
more information). Also shown as hash marks is the interannual standard deviation about the 
mean for each of the 19 products, b) monthly mean temperature for the New York City region, 
(°C), 1970-1999. The two observed products, the GCM ensemble average, and four points in the 
GCM distribution (lowest, 17th percentile, 83rd percentile, and highest) are shown. 

c) Mean annual precipitation for the New York City region, (cm), 1970-1999 in each of the 16 
GCMs, GCM ensemble, and Central Park observations. Also shown as hash marks is the 
interannual standard deviation about the mean for each of the 18 products, d) monthly mean 
precipitation for the New York City region, (cm), 1970-1999. Central park observations, the 
GCM ensemble average, and four points in the GCM distribution (lowest, 17th percentile, 83rd 
percentile, and highest) are shown. 
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1084 

1085 

1086 

1087 

1088 

1089 

1090 

1091 

1092 

1093 

1094 





Figure 4; Daily distribution (number of days per year) of: a) all-year mean, b) summer (June- 
August) maximum, and e) winter (Deeember-February) minimum temperature anomalies (“C), 
1980-1999 for Central Park observations (blaek line) and three GCMs (CSIRO, GISS, and MPI 
ECHAM5). 
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1095 

1096 

1097 

1098 

1099 

1100 

1101 

1102 

1103 

1104 

1105 

1106 

1107 

1108 

1109 

1110 

nil 

1112 

1113 


30 



Daily precipitation (mm) 


Figure 5; Daily distribution (number of days per year) of preeipitation (mm), 1980-1999 for 
Central Park observations (blaek line) and three GCMs (CSIRO, GISS, and MPI ECHAM5). The 
first bin, eontaining less than 10 mm, is not shown. 
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1114 


1115 

1116 

1117 

1118 

1119 

1120 
1121 
1122 

1123 

1124 

1125 

1126 

1127 

1128 



Year 



Year 


Figure 6; Combined observed (blaek line) and projected: a) temperature (°C) and b) annual 
precipitation (mm). Projected model changes through time are applied to the observed historical 
data. The three thick lines (red, green, and blue) show the ensemble average for each emissions 
scenario across the 16 GCMs. Shading shows the central 67 % range across the 16 GCMs and 3 
emissions scenarios. The bottom and top lines, respectively, show each year’s minimum and 
maximum projections across the suite of simulations. A ten-year filter has been applied to the 
observed data and model output. The dotted area between 2003 and 2015 represents the period 
that is not covered due to the smoothing procedure. 
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1129 

1130 

1131 

1132 

1133 

1134 

1135 

1136 

1137 

1138 

1139 

1140 

1141 

1142 

1143 

1144 

1145 

1146 

1147 


a) 



2020s 2050s 2080s 



2020s 2050s 2080s 


Figure 7: Seasonal a) temperature ehange (“C) and b) precipitation change (%) projections, 
relative to the 1970-1999 model baseline, based on 16 GCMs and 3 emissions scenarios. The 
maximum and minimum are shown as black horizontal lines; the central 67% of values are 
boxed, and the median is the thick line inside the boxes. 
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1148 


1149 



Daily winter minimum temperature (^C) 


1151 Figure 8; Daily distribution (number of days per year) of winter (Deeember-February) minimum 

1152 temperature anomalies (°C), for the New York Metropolitan Region in the CSIRO GCM. Blaek 

1 153 line, 1980-1999 hindeast; dotted line, 2080-2099 AIB seenario. 
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